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Abstract 

A formulation of Dirac's equation using complex-quaternionic co- 
ordinates appears to yield an enormous gain in formal elegance, as 
there is no longer any need to invoke Dirac matrices. This formula- 
tion, however, entails several peculiarities, which we investigate and 
attempt to interpret. 



1 Introduction 

One of the breakthroughs in the development of quantum mechanics was 
the discovery of a linear, relativistic wave equation for fermions by Dirac in 
1928, which is known as Dirac's equation P]I2] • This equation identified the 
electron spin as an intrinsic quantum number and predicted the existence 
of antiparticles. The corresponding Lagrangian is at the basis of quantum 
electrodynamics, and hence has played a major role on our path towards an 
understanding of elementary particles and their interactions. Dirac derived 
the equation, in a sense, by "taking the square root" of the Klein-Gordon 
equation, which is quadratic in time and space derivatives and the minimal 
version of a relativistically covariant wave equation. In order to overcome 
the technical obstacles in writing a linear equation, Dirac proposed that the 
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wave equation should be considered as a matrix equation for complex spinor- 
valued wave functions rather than complex scalar functions. The equation 
contains a set of 4 x 4 matrices, the so-called Dirac matrices, for which 
different representations satisfying the so-called Dirac algebra are possible. 

To the student who learns about these matrices for the first time, however, 
they may look somewhat arbitrary. One can easily derive the algebra and 
choose a specific representation, but nonetheless, on purely aesthetic grounds 
these matrices remain peculiar. A rather provocative statement by Einstein 
comes to mind: "When judging a physical theory, I ask myself whether I 
would have made the Universe in that way had I been God" . In the context 
of Dirac's equation, there is little sense in even questioning the correctness 
of the theory, but one may still ask whether there might not be a way to 
formulate it without availing oneself of these peculiar matrices. 

The structure of the matrices in Dirac's equation is intimately linked to 
relativistic covariance. Special relativity, however, can not only be formulated 
using the standard Lorentz four-vector notation, but also using a complex- 
quaternionic (CQ) parametrisation of space-time, as we have emphasised in a 
pedagogical review recently [3] . This parametrisation yields a significant gain 
in formal elegance, as Lorentz transformations are no longer implemented by 
multiplication with 4x4 matrices, but by multiplication with CQ numbers. 
Maxwell's four equations reduce to a single CQ equation. 

The question we wish to address in this article is whether one can for- 
mulate Dirac's equation (and quantum electrodynamics) in CQ coordinates, 
and to which extent such a formulation is equivalent to, or can be brought 
into a form equivalent to the standard theory. We will find that relativistic 
covariance implies immediately that the entities in the now two-component 
Dirac spinors are no longer complex numbers, but CQ numbers, and hence 
somewhat reminiscent of quaternionic quantum mechanics Among the 
peculiarities intrinsic to this formulation is an apparent doubling of solu- 
tions [30111, which we attribute to the existence of an additional, global 
SU{2) gauge symmetry, as well as the possibility to write solutions which 
carry spin one-half but do not possess a direction for the quantisation of 
this spin, which we eliminate from the spectrum of physical states through 
a suitable condition. 

Finally, we conclude that while it is indeed possible to gain an intriguing 
amount of elegance by formulating Dirac's equation in CQ coordinates, we 
have either failed to recognise a deeper physical principle at work or the price 
one has to pay for this gain in elegance is significant. 
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2 Special relativity 

We begin with a review of quaternions and the formulation of special rel- 
ativity within the framework of CQ numbers We introduce a complex 
algebra with generators 1, @ G C, such that 

@ 2 = -l (1) 

as well as a quaternionic algebra with generators k G H, such that 

= = Jfc* = -1, 

(2) 

r/ = — = k, jk = —kj = i, ki = —ik = j, 

which mutually commute: 

[ i ,@} = [ J ,@} = [k,@} = 0. (3) 

An arbitrary CQ number can be written as 

q = @t + ix + jy + kz 

with t, x,y,z G C. We further introduce a complex conjugate operation * 
which takes 

@ _> @* = _@ 

but leaves z, j, and unchanged, as well as a quaternionic conjugate operation 
~, which leaves @ unchanged but takes 

i — > i — —i, j —>■ j = — j, and k —>■ k = —k. 

Note that if Oi, o 2 G C ® H are two CQ numbers, the order of the product 
0\Oi is reversed under quaternionic conjugation only: 

(o x o 2 )* = o\o* 2 but o^=d 2 oi. (4) 

We further define a trace on the complex quaternions via 

tr : C ® H -> C, ? i-f - (g + q) , (5) 
which is C-linear and satisfies 

tr(pg) = tr(gp) (6) 
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for p, q G C <g> H. 

We label space-time (and other Lorentz contravariant quantities usually 
denoted by four vectors) by a purely imaginary CQ number 

q = @t + ix + jy + kz (7) 

with t,x,y,z G R, identify this subspace with Minkowski space, and de- 
note it by M. Complex conjugation * and quaternionic conjugation ~ then 
correspond to time reversal (T) and parity (P) transformations, respectively. 

The corresponding covariant quantity is given by its quaternionic conju- 
gate or parity reversed CQ number, 

q = @t — ix — jy — kz, (8) 

yielding the proper time interval 

j-2 2 2 2 /n\ 

—qq = —qq = t — x — y — z . (9 J 

Defining a scalar product on M 

(p,q) = ^(ip-q)(p-q)-(p + q)(p + q)) 

= -^{pQ + qp) 

= -^(pq + qp), (io) 

we find — qq = (q,q). Withp = @E+ip x +jp y + kp z and q = @t+ix+jy+kz, 
we have 

(p, q) = Et- p x x - p y y - p z z. (11) 

Note that @, k form an orthonormal basis of M. 

Let n = in x + jn y + kn z with nn = n x + n y + n 2 z = 1 and n x , n y , n z G R 
be a quaternionic imaginary unit vector. Then a Lorentz transformation is 
simply given by 

q - q' = uqu*, (12) 

where 

9 9 

u = e2 nd = cos - + n sin - (13) 
for a rotation by an angle 9 around n, or 

u = e^ n@A = cosh ^ + @n sinh ^ (14) 



for a boost by a Lorentz angle A in direction n. 

Clearly the covariant CQ number q transforms as 

q — > q' = u*quj. (15) 

With ijju = ujuj = 1, the Lorentz invariance of the scalar product (p, q) is 
evident. 



3 Dirac equation 

We now apply the complex-quaternionic formulation of special relativity to 
relativistic wave equations. We begin with the simplest CclSG, 1.6. ? db free scalar 
field. Let <&(q) be a function of one purely imaginary CQ number q G M, 
which transforms as a scalar under Lorentz transformations: 

$(<Z)^$V) = *(?)• (16) 
The contravariant differentiation operator is defined by 

d = @^ -i^L - j^L - = m t -id x -jd y - kd z (17) 

and transforms according to 

D —> D' = luDlj*. (18) 

Note that —Dq = 4, which is equivalent to d^x^ = 4 in standard notation. 
With the Lorentz invariant operator —DD the free Klein-Gordon equation 
reads 

{—DD + m 2 )&(q) = 0, (19) 

where the real scalar m represents the particle mass. The generally known 
solutions of (|19|) may be written as 

$( g )= e ®<^>, (20) 

where the energy-momentum p = @E + ip x + jp y + kp z satisfies 

-pp = E 2 - p 2 x - p\ - p\ = m 2 . (21) 
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In principle, one could also write down solutions of the form 



$( g ) = e nM , 

where n is an arbitrary quaternionic imaginary unit vector. These solutions, 
however, do not correspond to any physically new states and, in contrast to 
(J20J), are not ordinary complex functions. We hence give them no further 
consideration. 

Within this framework, Dirac's equation for the two-component field 
Oi,^2) T is given by 



-m D \ / ifti 
D -m J \ i>2 



0. (22) 



This equation is covariant if m is a real scalar and the component fields 
transform like 

ip[ = uifji and ip' 2 = u*if)2. (23) 

This shows that ipi^ have to be CQ-valued in general. Therefore, Lorentz 
transformations are represented by diagonal matrices (|23|) . which is clearly a 
simplification compared to the usual spinor representations used in Dirac's 
theory [2|. 

Iteration of (}2*2*|) yields the Klein-Gordon equation (|TT?|) for every compo- 
nent ipi 2 and therefore the relativistic dispersion 



E' z = p 2 + m 2 , 



where p is the spatial momentum. From (|23|) for rotations one finds that, as 
u(2rc) = — 1, particles obeying (J22j) have spin one-half. 

Solving Dirac's equation in the rest frame (p = 0) we obtain the solutions 

* = ( ] e T @mt , ijj G C ® H, (24) 



where the upper sign holds for particles (E > 0) and the lower sign for 
antiparticles (E < 0). If one regards C <8> HI as a four- dimensional complex 
vector space, (|2^j) represents eight linearly independent solutions of (j22J), i.e., 
twice the number as for the ordinary Dirac equation 2j. This doubling of 
solutions has led Edmonds |5 a and Gough jS] to investigate the possibility 
of the existence of hidden quantum numbers and associated observables. Dc 
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Leo |7J proposed that one may reduce the number of solutions by transform- 
ing ([22)1 into a one-component equation which then requires left and right 
multiplications by quaternions. We show in Section 6 below that there ex- 
ists an additional symmetry of ([22)1 . which reduces the number of linearly 
independent solutions from eight to four through gauge invariance. 

Solutions with finite velocity are obtained by writing the exponent in (|24j) 
Lorentz scalar. The Ansatz 

*± = ( ) (25) 



yields the linear system of equations 



m ±@p \ / ip±(p) 
±@p m J \ £±(p) 



(26) 



This implies the spinor condition 



e±(p) = T~@ P MP) (27) 
m 

and the energy-momentum relation 

pp + m 2 = -E 2 +p 2 + m 2 = 0. (28) 
Note that (|27j) is invariant under Lorentz transformations and is satisfied by 



4 Spin and spin directions 

In the previous section we have seen that the solutions of ()22|) describe par- 
ticles with spin one-half. For further investigations note that, as the spin 
operator is the generator of rotations in the rest frame, we have to restrict 
ourselves to particles in their rest frame, which are represented by states of 
the form (|24|) . Let us first identify the components of the spin operator. 
Consider a rotation around the x-axis by an angle 9. With uj = uj* we obtain 

W>i,2 = = e- 6 *"^- (29) 

The components of the spin operator are now easily read off to be 

@ @ @ 

S x = S v = TTJj and S z = —k. (30) 
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With = 2k and cyclic permutations one readily verifies that (J3U|) satisfy 
the angular momentum algebra 



[S X , Sy] = @S Z , [Sy, S Z ] = @S X , [S z , S X ] = @Sy. (31) 

For the square of the spin operator we obtain S 2 = |, as expected. The spin 
operator S is hence implemented through a simple left multiplication with 
the quaternionic basis elements. Our construction differs from that of De 
Leo |2] and Gough 9j, who proposed operators consisting of simultaneous 
left and right multiplication with k. 

Let us now choose the z-axis as spin quantisation direction. The solutions 
of (J22J) which are simultaneously eigenstates of S z with eigenvalues m z — ±g 
are given by 

m z = \: (1 + @fc)* , (i + @j)^o 

m z = -\ : (@i + j)*o, (~@-k)%, ( ° _) 



where 

' 1 
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for particles. Antiparticle solutions are constructed similarly. With the usual 
definition for the raising and lowering operators = S% ± ©5*2 one easily 
verifies that the two subspaces 

span c {(l + @fc)^ , (@i+j)ty } and span c {(i+@j)* , (_@_fc)* } (33) 

are closed under the spin algebra (|31|). Note that we can convert the spin 
eigenvectors (}3*2j) for the same eigenvalue into each other by right multipli- 
cation with quaternionic basis elements, e.g., 

{l + @k)ty i = {i + @j)^ and (@i + j)V i = (-@ - fc)tf„. (34) 

The right multiplication with the other basis elements j, k does not yield 
independent states. Hence we can connect the two spin eigenspaces (JUnj) by 
right multiplication with i,j, k. We will return to this issue in Section 6. 

Another peculiarity of the CQ formulation of Dirac's equation is the fol- 
lowing. The states 

tf , i*o, i*o, k^ (35) 
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form a basis of the subspace of particle solutions. They are somewhat pe- 
culiar, however, as they mix the two spin eigenspaces (J23), e.g., we have 



^o = ^(l + @A;)tto+®(-@-fc)*o- (36) 

Note that the eigenvectors (f3*2*J) of S z satisfy the condition \I/ T \I/ = (where 
T denotes transposition) whereas the states (|33|) do not. In general, the 
following two statements are equivalent: 

^ T il> = o (g. S)V = ±i*, (37) 

where e is a unit vector in R 3 . The dot • denotes the standard scalar product. 
The right hand side of ()37|) states that the spin of the particle points in the 
direction e. The states are solutions to Dirac's equation and possess 
spin one-half, but in contrast to the usual situation there is no direction e of 
the spin. Rewriting (|37jl for the states ({231 we obtain the condition 

^ = (e • S^if) = ±-ip. (38) 

With ifj = vpo + iifii + jifj 2 + kips, the equation on the left can be rewritten 

^ + ^ + ^2+^3 = 0. (39) 

Since C®H is isomorphic to the space of complex 2x2 matrices |10j . (}3*9"j) is 
equivalent to the statement that ij) is not invertible in C ® HI. The condition 
(jBHt holds for antiparticles as well. 

To prove (J37|) we start with the eqivalent condition (J38)) . First note that 
e ■ S = |n with n G H\R. Then multiplying ()38|) by tJ) from the right and 
assuming ipip ^ 0, we obtain the contradiction 

1 @ 

2 2^ 

Hence, in order to fulfil ()38|) we need ipif) = 0. 

The other direction is proven by rewriting ip G C <S> EI as 

ip = a + @b + cm + @dm' (40) 

with m, m' quaternionic imaginary unit vectors and a, b,c,d G K. Then, by 
setting 6 = 0, we have the condition 

^■0 = (a + cm) (a — cm) — (a + cm)@dm' + @dm(a — cm) — d 2 = 0, 
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which yields 

a 2 + c 2 - d 2 = 0, (mm' + m'm)cd = 0. (41) 

For ip ^ the equation on the left implies d ^ 0. In the first case we 
have c = and can choose n = m! in ()38|) . in the second case we have 
mm' + m'm = and we can choose n = ±^(a + cm)m', which completes the 
proof. 

The above implies that states obeying (|H7j) may be interpreted as ordinary 
spin one-half particles with usual spin properties. There are, however, two 
subspaces ()33j) closed under the spin algebra ()31|) rather than one, which can 
be connected by right multiplication with the quaternionic basis elements. In 
addition to these familiar states there are also states with ^ T ^> ^ which are 
not interpretable as ordinary particles. We will return to these peculiarities 
in Section 6. 

5 Quantum electrodynamics 

So far we considered only free particles. We now turn to Dirac fermions 
coupled to an external electromagnetic potential, which we describe with 
the Lagrangian formalism. For free Dirac fermions, the Lorentz invariant 
Lagrangian density is given by 

^ = *((*'*)("S -m)(t)) 

= tr + ^ 2 *£>^2 - m$Va - W^i) , (42) 

where the trace is as defined in JSJ). The free action is given as usual by the 
integral 

S = [ C d\ (43) 
Jm 

where d\ denotes the volume element in M = IR 4 . With D* = —D, D* = —D 
we find 

£* = tr (-^iltyl - $2Dip* - m$iV>2 - m^ipf) , (44) 

which with integration by parts yields the reality condition <Sq = «So- The 
trace in (}4~2|) implies S = S such that Sq G R. 
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To couple the Dirac fields to the electromagnetic potential we require 
invariance under the local gauge transformation 

_> = ^ e -@e«(9) (45) 

with e a real constant and a:M->la scalar function. This implies the 
minimal coupling procedure D — > D + @eA in (|22j) and the transformation 
rules for the vector field A = + %A X + jA y + kA z : 

A — > uAuj* under Lorentz transformations, (46) 
A — > A + Da under £7(1) gauge transformations. (47) 

The interaction is hence given by the Lagrangian density 

£ int = tr (@e^M + @e^^ 2 ) , (48) 

which is real as ^4* = —A and A* = —A. We identify e with the electric 
charge. The electromagnetic field strength for A is given by |3] 

F=~(dA-DA), (49) 

and the Lagrangian density for the free electromagnetic field is 

C A = \{F 2 + (F*) 2 ) , (50) 

which is real since F = —F. The Lagrangian density of quantum electrody- 
namics is hence 

£QED = £o + Ant + f~-A- (51) 

From the non-relativistic limit of Dirac's equation, we can further read off 
the g-f actor of the electron to equal g — 2. 

6 Quaternionic gauge invariance and the con- 
nection to the conventional Dirac equation 

We have seen in Section 3 that a doubling of solutions occurs in the CQ 
version of Dirac's equation. This is also evident from the existence of two 
closed spin eigenspaces f!33f) rather than one. We have seen that these two 
eigenspaces are connected by right multiplication with k. 
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We generalize these transformations to 

$ _> & = m e - nP , (52) 

where n = in x + jn y + fcn 2 is a quaternionic imaginary unit vector and 
/3 6 M. Note that e _ri/3 can be represented in the basis 1, i, j, k of H (see also 
(HHJ)). The Lagrangian density (J4*2~j) is, by use of invariant under this 
transformation. This is, however, not a physical symmetry of the system 
but rather an invariance of description within the CQ formulation of Dirac's 
equation, as we will explain now. 

Note first that the unit quaternion e~ n/3 represents a point on the three 
sphere S 3 C 1R 4 . Therefore, ()52|) represents a global SU(2) gauge symme- 
try. The standard procedure to introduce a gauge symmetry is to require 
the system to be invariant under a local gauge transformation. To meet this 
requirement one introduces a gauge field (like A for the U{1) gauge trans- 
formations in Section 5) which is minimally coupled to the matter field 
This coupling represents the interaction of the matter with the gauge field. 
The extension of (|52J) to a local gauge transformation, i.e., (3 = /3(q), how- 
ever, does not seem possible within the CQ formulation of Dirac's equation 
presented here. To see this, consider for simplicity (3{q) = x and introduce 
a gauge field B by replacing the differentiation operator D by D + B. Then 
we have to require 

((D + B)V)' = (D + B)'^' = D^' + BW (53) 

under the (local) transformation (p)2*j) . For j3{q) = x this yields 

(Dm)e~ nx + B^e~ nx = (D^)e- nx + i^ne~ nx + B'^e~ nx , (54) 

where we have used the definition of the differentiation operator D (see (|17|)). 
As n G H and there is no general rule to express \l/n in terms of n\I/, it is 
impossible to satisfy (p>fj) with B a CQ- valued function or matrix jllj . 

It is hence not possible to elevate the global "gauge" symmetry (jo^j) to a 
true and local gauge symmetry. The symmetry is, however, a gauge symme- 
try, as it corresponds to an invariance of our description (as gauge symmetries 
generally do) rather than an invariance of the physical system. It is signif- 
icant for the interpretation of the CQ formulation of Dirac's equation, as it 
connects the two spin eigenspaces (J3l!j) but leaves the spin eigenvalues invari- 
ant. Identifying states connected by quaternionic gauge transformations (|52j) 
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with each other, the number of particle solutions of the free Dirac equation 
(J22J) reduces from two to one for each spin eigenvalue. Note that (J52|) is 
trivially extended to quantum electrodynamics with A —>■ A. 

Recalling (|3*5jl. we see that states with ^ T ^> ^ correspond to super- 
positions of states formulated in different quaternionic gauges. The crucial 
point is that no gauge transformation exists which transforms a state with 
ty T ty 7^ into a state which belongs to only one of the subspaces (J3l?|) . e.g., 
there is no solution for n and (3 satisfying 

# oe r n £ = a (l + @k)V + b(@i + i)* , a,beC. (55) 

By contrast, spin eigenstates can be transformed by (|52|) to states which 
belong to any one of the spin eigenspaces (J3l?|) . For example, 

-=(1 + @k)V + -^=(i + @j)*o = (1 + @fc)* e^. (56) 
v 2 v 2 

This does, however, not provide us with a definite answer to the question 
whether states with \I/ Tl ^ ^ should be interpreted physically or should be 
excluded from the theory by requiring ^ T ^> = for physical states explicitly. 
The contemplations summarized above suggest the latter. 

To verify that there are no further gauge symmetries, note that the La- 
grangian density (|42jl is invariant under right multiplication of a CQ number 
E only if EE* = 1. Then the statement 

SS* = 1 =>■ E = eg, (57) 

where c £ C, g G H, and |c||g| = 1, shows that there are no gauge symmetries 
besides the U(l) symmetry as described in Section 5 and the quaternionic 
gauge transformations (|52|). The proof of (|57|) is straightforward, but some- 
what to technical for this article. 

We conclude this section by relating the CQ version of Dirac's equation 
(J 12 12 j) to the standard formulation with complex matrices j2]. For this purpose, 
we expand the CQ spinors ipip in up- and down-spin components, where we 
can restrict ourselves due to (|52~j) (at least for states satisfying ty T ty = 0) to 
the first of the two subspaces (HUD), i.e., we write 

ipx = (l + @fc)ci + + j)c 2 , ip 2 = (1 + @k)c 3 + (@i + j)c 4 (58) 

with d : M — > C. Inserting this in (|2*2^) yields four independent equations 
which combine to 

(<&fd li -m)C = with C = ( Cl ,c 2 ,C3,c 4 ) T , (59) 
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where the gamma matrices are given in the chiral representation: 

with and 1 are the 2x2 zero and unit matrix and a 1 the Pauli matrices. 
The CQ formulation with \E fTl $ r = presented here is hence equivalent to the 
standard complex formulation of Dirac's equation. 

7 Conclusion 

In this article, we have extended the complex-quaternionic formulation of 
special relativity [3] to Dirac's equation, obtaining an equation for a two- 
component CQ field. In this formulation there is no need to invoke Dirac 
matrices, and Lorentz transformations are given by diagonal matrices rather 
than complicated spinor representations. The price to pay for these sim- 
plifications is, however, outrageous. First, we encountered a doubling of 
solutions, a problem which we overcame by attributing them to different 
gauge choices of a global 577(2) gauge symmetry. Second, we encountered 
particle solutions with spin one-half but without a spin quantisation direc- 
tion. According to conventional quantum mechanics these solutions appear 
to be unphysical, but this did not emerge from the formalism we developed 
above. This led us to derive a condition to ban them from the Hilbert space, 
which we imposed a posteriori. Furthermore, there are mathematical prob- 
lems inherent to the formulation due to the fact that C ® H is not a division 
algebra 0JH2|. 

In summary, we found no indication that the CQ formulation of Dirac's 
equation would lead to a deepening of our understanding of relativistic fer- 
mions. In our opinion, the problems we encountered with this formulation 
eradicate the gain in formal elegance we were hoping to obtain. 
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